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Abstract 


Background: Human coronaviruses (HCoVs) OC43, 229E, NL63 and HKU1 cause upper and 
lower respiratory tract infections. HCoVs are globally distributed and the predominant species may 
vary by region or year. Prior studies have shown seasonal patterns of HCoV species and annual 
variation in species prevalence but national circulation patterns in the US have not yet been 


described. 


Objectives: To describe circulation patterns of HCoV species OC43, 229E, NL63 and HKU1 in 


the US. 


Study Design: We reviewed real-time reverse transcription polymerase chain reaction (rRT-PCR) 
tests for HCoV-OC43, -229E, -NL63 and -HKU1 reported to The National Respiratory and Enteric 
Virus Surveillance System (NREVSS) by U.S. laboratories from July, 2014—June, 2017. We 
calculated the total number of tests, and percent positive by week. For a subset of HCoV positive 
submissions with age and sex of the patient available we tested for differences in age and sex 
across the four HCoV_ species using Chi Square and Kruskal Wallace tests. 
Results: 117 laboratories reported 854,575 HCoV tests; 2.2% were positive for HCoV-OC43, 
1.0% for HCoV-NL63, 0.8% for HCoV-229E, and 0.6% for HCoV-HKU1. The percentage of 
positive tests peaked during December — March each year. No significant differences in sex were 
seen across species, although a significant difference in age distribution was noted. 
Conclusions: Common HCoVs may have annual peaks of circulation in winter months in the US 
and individual coronaviruses may show variable circulation from year to year. Different HCoV 
species may be detected more frequently in different age groups. Further years of data are needed 


to better understand patterns of activity for HCoVs. 
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Background 


Human coronaviruses (HCoV) HCoV-NL63, HCoV-HKU1, HCoV-229E, and HCoV- 
OC43 circulate worldwide and cause a range of respiratory symptoms [1]. Infections are often 
asymptomatic or associated with mild to moderate upper respiratory tract illness in 
immunocompetent children and adults; HCoVs are considered the second most common cause of 
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the common cold [2]. Infections can also result in lower respiratory tract illness including 
bronchiolitis and pneumonia, especially in immunocompromised individuals, infants, and older 
adults [1]. Increased availability of molecular test methods and more frequent testing for multiple 
respiratory pathogens have allowed for opportunities to characterize circulation patterns of 


individual HCoVs. 


Although HCoVs are globally distributed, the predominant species may vary by region or 
year [3,4,5]. Previous studies have shown seasonal patterns of HCoV species and annual variation 
in species prevalence [3,4,6]. However, national circulation patterns across the United States have 
not been described and few studies have described circulation of all four HCoVs across multiple 


years [6]. 


Objectives 


Our objective was to describe laboratory detections of HCoVs -NL63, -HKU1, -229E, and -OC43 
in the United States during 2014-2017, using data collected by The National Respiratory and 


Enteric Virus Surveillance System (NREVSS) [7]. 


Study Design 


NREVSS is a passive surveillance network established by the U.S. Centers for Disease 
Control and Prevention (CDC) in the 1980s that collects specimen test results for several 
respiratory and enteric viruses from multiple laboratories across the United States [7]. NREVSS 
currently collects data from three different sources: directly from clinical, state, and local 
laboratories; indirectly from state or local partners on behalf of laboratories within their 
jurisdictions; and indirectly through the Public Health Laboratory Interoperability Project 


(PHLIP). 


PHLIP is a collaborative partnership between CDC, state and local public health labs 
(PHLs), and the Association of Public Health Laboratories (APHL) to strengthen the submission 
of automated specimen-level surveillance laboratory test results directly to CDC [8]. Laboratories 
submitting to NREVSS via PHLIP submit specimen level results for respiratory virus tests along 
with patient demographics, such as sex and age. The remaining majority of NREVSS participants 
(i.e. non-PHLIP submitters) submit weekly aggregates of positive detections for the four HCoV 


species by RT-PCR along with the aggregate number of RT-PCR HCoV tests performed. 


To better understand HCoV circulation in the US, we assessed reports of specimens tested 
for HCoVs, and submitted to NREVSS during July 1, 2014 — June 30, 2017. Reporting laboratories 
were excluded if they did not report HCoV results at the species level. We summarized the total 
number of HCoV tests submitted during the study period and calculated the overall percent 
positive for each HCoV species by week. Test numbers and positive HCoV results were also 


summarized separately for individual US census regions. 


To understand demographics of patients with HCoV detections we further analyzed the 
subset of NREVSS data reported through PHLIP. We selected reports of specimens tested for all 
four HCoV species and calculated the percent positive for each HCoV species. We used all PHLIP 
reports with a single positive HCoV detection to test differences in age distribution among the four 
HCoV species using the Kruskal Wallis Test. Differences in sex distribution among the four HCoV 
species were tested using the Chi Square Test. We then summarized viral co-detections for each 
HCoV species using reports of specimens tested for all four HCoVs and which were also tested 
for the following viruses: parainfluenza viruses 1-4, respiratory syncytial virus (RSV), human 
metapneumovirus, human adenovirus, rhinovirus/enterovirus, influenza A and influenza B. 


Analysis was performed using R version 3.3.1. 


Results 


During July 1, 2014— June 30, 2017, 854,575 HCoV tests were reported by 117 laboratories 
in 42 states submitting to NREVSS. Overall, 18,804 (2.2%) were positive for HCoV-OC43, 8,558 
(1.0%) for HCoV-NL63, 7,001 (0.8%) for HCoV-229E and 5,225 (0.6%) for HCoV-HKU1. The 
number of HCoV tests submitted to NREVSS per week varied seasonally, with the testing peak 
occurring each year in winter, generally between December and March (Figure 1A). Overall HCoV 
testing increased during the three years (Figure 1A). The percent of HCoV positive tests varied 
throughout each year, and also peaked each year between December and March. The percent 
positive varied annually by HCoV species (Figure 1B). HCoV-OC43 demonstrated a distinct peak 
each of the three years, with a less pronounced peak in 2016. HCoV-NL63 and HCoV-HKU1 
demonstrated similar patterns to one another; both had a small peak in 2015 and larger peaks in 
2016, although only HCoV-NL63 had a small peak in 2017. HCoV-229E showed a slight peak in 
2015, no peak in 2016 and a relatively large peak in 2017 (Figure 1B). The highest percent positive 
for any single species was 7.6% of tests positive for HCoV-OC43 in the week beginning December 
31, 2016. Across each census region, minimal differences in seasonal and annual patterns of 
percentage of tests positive for each HCoV species were seen compared to national data (Figure 
2). The most notable difference was in the percentage of positive HCoV-OC43 tests during the 
2016-2017 season, with the West region showing a peak percent positive of 5.0% and the Midwest 


region showing a peak percent positive of 12.4%. 


Data reported to NREVSS through PHLIP was further analyzed to understand sex and age 
characteristics. During the study period, 20,806 specimens tested for all four HCoVs were 
submitted via PHLIP from six laboratories. Overall 1,569 tests (7.5%) were positive for any 
HCoV; 852 (4.1%) were positive for HCoV-OC43, 255 (1.2%) for HCoV-NL63, 335 (1.6%) for 
HCoV-229E and 154 (0.7%) for HCoV-HKU1. The majority of specimens with a single HCoV 
detection (92.2% of 1543 specimens) included the sex of the patient, and approximately half 
(50.6%) of all these HCoV detections were reported as male. No significant difference was seen 
in the sex distribution between the four HCoVs (p=0.19). Age was available for 1,016 (67%) of 
specimens with a single HCoV detection and the median patient age of these specimens was 23 
years (range 0-96 years). The patient age distribution of specimens with a single HCoV detection 
was significantly different between species (p<0.01) (Figure 3). The median ages of patients with 
specimens testing positive for a single HCoV species were as follows: HCoV-OC43, 24 years; 
HCoV-NL63, 11 years; HCoV-229E, 30 years; and HCoV-HKU1, 19 years. For HCoV-OC43, 
HCoV-NL63, and HCoV-HKU1 >45% of detections were in children <18 years old (Figure 3). 


By contrast, 31% of HCoV-229E detections were in children <18 years old (Figure 3). 


Among the 1,569 HCoV positive detections reported via PHLIP, 1,538 (98%) were also 
tested for parainfluenza viruses 1-4, respiratory syncytial virus (RSV), human metapneumovirus, 
human adenovirus, rhinovirus/enterovirus, influenza A and influenza B. Among these, 68.6% 
reported a single HCoV species detection only, 1.7% reported two or more HCoV species, and 
30.2% detected another respiratory virus. The most common HCoV co-detections were HCoV- 
OC43 with HCoV-NL63 (8 specimens, 0.5%), and HCoV-OC43 with HCoV-229E (8 specimens, 
0.5%). The most common co-detected non-HCoV viruses were RSV (11% of HCoV positive 


specimens), rhinovirus/enterovirus (6.6%), and influenza A (5.7%); 51 (3.3%) specimens had >2 


viral species detected in addition to HCoV. Co-detection patterns were broadly similar among the 


four HCoVs (Table 1). 


Discussion 


This report is the first to describe the national patterns of circulation of the four common 
HCoV species in the United States during a multi-year period. During the study period, HCoVs 
showed a peak prevalence during December— March each year, which coincides with the winter 
respiratory virus season [7,9]. HCoV-OC43 was the most commonly detected HCoV with 2.2% 
of all tests positive. Different HCoV species predominated in different years; HCoV-OC43 
appeared to peak annually, while HCoV NL-63, HCoV-HKU1, and HCoV-229E showed more 
variability, with distinct peaks in one or two of the three years studied. This is consistent with 
previously published site-specific data indicating that individual species may only demonstrate 


peak activity every 2-3 years [10,11]. 


Factors associated with annual differences in activity for HCoV species are currently 
unknown. Individual HCoV species activity could fluctuate independently, or cross-immunity 
within or between Alphacoronavirinae (HCoV-229E and HCoV-NL63) and Betacoronavirinae 
(HCoV-HKU1 and HCoV-OC43) might affect annual activity of the four HCoVs [4,12]. During 


the study period, HCoV-229E and HCoV-NL63 did not show large contemporaneous peaks of 


activity, although both showed smaller peaks of activity in 2014-2015. When data was visualized 


by census region, annual and seasonal patterns were similar to those seen nationally (Figure 2). 


The age distribution of patients with reported HCoV infections differed between HCoV 
species (Figure 3). HCoV-229E detections were more common in adults >18 years old compared 
to HCoV-HKU1, -NL63, and -OC43. HCoV-229E has previously been reported as 
disproportionately affecting immunocompromised individuals relative to the other HCoV species 
[4], possibly affecting the median reported age at infection. HCoV-NL63 showed the lowest 
median age of infection, and has been shown to be associated with croup in young children [13], 
which may lower the average age of infection. Non-HCoV viral co-detections were seen in 30% 
of specimens positive for HCoV, and 3.3% of specimens had two or more co-detected viral species. 
The clinical impact of coronaviruses in co-detections is not fully understood, with prior studies 
reporting both increased and unchanged morbidity and mortality with respiratory viral co- 
detections [14]. Single infections with HCoVs have been associated with morbidity due to lower 
and upper respiratory tract infections [4], however additional studies have indicated that HCoV 
infections in children may not be associated with hospitalization or increased severity of illness 


[15]. 


There were limitations to this report. NREVSS is a passive, voluntary surveillance system, 
collecting results from specimens submitted to U.S. laboratories. Many HCoV infections are 
subclinical or mild, and do not require clinical care; therefore these infections are unlikely to 
require laboratory testing and would not be captured by NREVSS. The relative proportions of 
HCoV species reported here may not be representative of all HCoV infections. Within NREVSS 
all reporters of HCoV surveillance data at the species level were included, including those 


inconsistently reporting over time. This may result in certain laboratories or regions being 


overrepresented at certain times e.g. during winter, when the majority of respiratory virus tests are 


conducted. 


Data reported to NREVSS through PHLIP represent a smaller subset of six laboratories 
that also report additional data including sex and age; data reported through PHLIP may not be 
representative of the entire NREVSS population. The overall percent of HCoV positive tests varied 
between PHLIP and NREVSS; the percent of HCoV positive tests submitted to NREVSS via 
PHLIP was 7.5%, vs 4.6% for NREVSS as a whole. This may result from qualitative differences 
in reporting laboratories; reporting through PHLIP is limited to military, state and public health 
laboratories, whereas data submitted directly to NREVSS through non-PHLIP sources primarily 
includes clinical laboratories. State and local public health laboratories test for a variety of reasons 
including surveillance and public health response to outbreaks and clusters of cases, whereas 


clinical laboratories often test in the course of managing individual cases. 


Within the peak in detections noted during December-March, we were not able to define 
precise seasonal onset and offset periods for any individual HCoV species. Additionally, we 
assessed HCoV percent positivity over time based on the date of submission of a report to 
NREVSS rather than the date of specimen collection, and we anticipate that our findings may 
reflect this slight delay. Finally, aggregate data reported to NREVSS via PHLIP might include 
multiple specimens from the same patient, potentially impacting the demographic characteristics 


of reported cases. 


Surveillance of HCoVs is important to determine seasonality and annual circulation 
patterns. Continued use of respiratory virus multiplex assay panels could facilitate further 


definition of HCoV circulation through public health surveillance in the future. Further years of 
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data are needed to better understand patterns of activity for HCoVs, and studies with additional 


epidemiologic data will be useful to better characterize HCoV burden and spectrum of illness. 


Disclaimer: The findings and conclusions in this report are those of the authors and do not 


necessarily represent the official position of the U.S. Centers for Disease Control and Prevention. 
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Figure 1: A) The number of tests performed to detect any of the four HCoVs OC43, NL63, 229E, 
and HKU1 reported to the National Respiratory and Enteric Viruses Surveillance System 
(NREVSS) by week, July 2014 — June 2017, B) The percentage of tests positive for HCoVs OC43, 


NL63, 229E, and HKU1 reported to NREVSS by week, July 2014 — June 2017. 
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Figure 2: The percentage of tests positive for coronaviruses OC43, NL63, 229E, and HKU1 
reported to the National Respiratory and Enteric Viruses Surveillance System (NREVSS) by week, 


stratified by US Census Region, July 2014 — June 2017. 
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Figure 3: Age distribution of HCoV- 229E, HKU1, NL63 and OC43 positive tests with age 
available reported to NREVSS via the Public Health Laboratory Interoperability Project (PHLIP), 
from July 2014 — June 2017. Note, specimens with more than one coronavirus detected were 


excluded (n=18). 
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Table 1: Frequency and percentage (in parentheses) of viral co-detections for by individual 
HCoV detection and all HCoV detections, reported to NREVSS via PHLIP from July 2014 — 























June 2017 
> 
© 
= 
(oe) (oe) = L 
S 8 R or sE 
Se. Qro m s z 3 
S § § ag i BSE of 
mom ea n Æ E E 48 
HCoV- 836 8 8 4 28 100 23 54 59 44 10 259 
OC43 n= (100) (1.0) (1.0) (0.5) (3.3) (12.0 (2.8) (65) (7.1) (5.3) (1.2) (31.0 
836 ) ) 
HCovV- 8 253 3 5 5 26 10 13 19 11 0 74 
NL63 (3.2) (100) (1.2) (2.0) (2.0) (0.3 (4.0) (5.1) (7.5) (43) (0.0) (29.2 
n=253 ) ) 
HCoV- 8 3 325 0 16 21 9 11 11 27 5 89 
229E (2.5) (0.9) (100) (0.0) (4.9) (65) (2.8) (3.4) (4) (8.3) (1.5) (27.4 
n=325 ) 
HCovV- 4 5 0 151 5 21 5 10 13 8 2 53 
HKU1 (2.6) (3.3) (0.0) (100) (3.3) (13.9 (3.3) (6.6) (86) (5.3) (1.3) (35.1 
n=151 ) ) 
All 836 253 325 151 54 164 47 86 101 87 17 465 
HCoVs (544 (164 (21.1 (9.8) (3.5) (0.7 (8.1) (5.6) (66) (5.7) (1.1) (0.2 
) ) ) ) ) 
= respiratory syncytial virus, HMPV = human 


* PIV = parainfluenza virus, RSV 


metapneumovirus, HAdV = human adenovirus, RV/EV = rhinovirus/enterovirus, Flu A = 


influenza A, Flu B = influenza B 
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